Isotropic contact forces in arbitrary representation: 
Heterogeneous few-body problems and low dimensions 
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The Bethe-Peierls asymptotic approach which models pairwise short-range forces by contact con- 
ditions is introduced in arbitrary representation for spatial dimensions less than or equal to 3. The 
formalism is applied in various situations and emphasis is put on the momentum representation. In 
the presence of a transverse harmonic confinement, dimensional reduction toward two-dimensional 
(2D) or one-dimensional (ID) physics is derived within this formalism. The energy theorem relating 
the mean energy of an interacting system to the asymptotic behavior of the one-particle density 
matrix illustrates the method in its second quantized form. Integral equations that encapsulate the 
Bethe-Peierls contact condition for few-body systems are derived. In three dimensions, for three- 
body systems supporting Efimov states, a nodal condition is introduced in order to obtain universal 
results from the Skorniakov Ter-Martirosian equation and the Thomas collapse is avoided. Four- 
body bound state eigenequations are derived and the 2D '3-1-1' bosonic ground state is computed 
as a function of the mass ratio. 

PACS numbers: 03.65.Nk,03.65.Ge,21.45.-v,34.50.-s 



I. INTRODUCTION 

The non-perturbative zero-range limit of a pairwise in- 
teraction was introduced in the context of nuclear physics 
for the description of the deuteron [J-Q- I^ ^^^^ mod- 
eling, the finite range pairwise interacting potential is 
replaced by a zero-range potential and a peculiar asymp- 
totic behavior is imposed on the wave function as the 
two interacting particles approach each other. More than 
60 years after these pioneering works, the Bethe-Peierls 
approach appears to be especially relevant for modeling 
pairwise interaction in few- and many-body systems of 
ultra-cold atoms where details of the interaction are use- 
less for understanding a large class of low-energy pro- 
cesses |4l-[l7|. This remarkable simplicity follows essen- 
tially from two features: first, the temperature in these 
systems is sufficiently low that s-wave scattering is the 
dominant process in two-body collisions [185; second, the 
three-dimensional (3D) s wave scattering length (denoted 
as) in the two-body scattering can be tuned by use of a 
magnetic Feshbach resonance (FR) toward an arbitrary 
large value with respect to the range of the interatomic 
forces (denoted by 63) and which is of the order of the 
van der Waals radius [l^, 1201 : 
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where Ce is the van der Waals coefficient of the inter- 
atomic potential and /z is the reduced mass of the two 
interacting particles. In the vicinity of a broad FR the 
two-body scattering cross section Ana^ is very large with 
respect to the non-resonant cross section (^ 47163) ''^^'^ 
this justifies the use of the 3D Bethe-Peierls approach 
only parameterized by 03. 

Ultra-cold atoms in reduced geometries are the sub- 
ject of intensive experimental and theoretical studies [21j . 



Decisive progress have been obtained both in quasi-lD 
systems (for exam ple with the observation of the Tonks 
Girardeau gas [22l, 123]) and in quasi-2D systems. A gen- 
eral review of low-dimensional systems can be found for 
example in Refs. [24s, 25]. Systems in reduced dimensions 
are achieved by using very anisotropic trapping poten- 
tials leading to a _D-dimensional behavior where D = 1 
(quasi-lD systems) 01 D = 2 (quasi-2D systems). In 
an ideal D-dimensional atomic waveguide, noninter act- 
ing atoms move freely along D direction(s) while they 
remain frozen in the lowest state of the transverse zero- 
point motion characterized by a length denoted as atrans- 
The associated energy E'trans = f^ I {'^'IJ-a'tra.ns) defines the 
limit of energy beyond which higher transverse states are 
populated and a 3D-like behavior is progressively recov- 
ered for increasing energies. In the actual experimen- 
tal state-of-the-art, atomic waveguides have a transverse 
atomic length which is large as compared to the 3D po- 
tential radius (63). Therefore for low energy processes 
i.e., for collisional energies E such that \E\ <C fi.^/(/i6|), 
the short-range pairwise interacting potential can be de- 
scribed through the 3D Bethe-Peierls approach. For col- 
lisional energies much smaller than -Etrans the transverse 
excited states of the trap are populated only via vir- 
tual processes and for large interatomic separations r 
(r 3> fltrans) Only the ground transverse state is occupied. 
In this regime of collisional energies and in the free D- 
dimensional subspace, scattering process can be deduced 
from a D-dimensional effective pairwise interaction char- 
acterized by a finite range denoted bjj, which is of the 
order of the transverse length atrans- Hence for low en- 
ergy D-dimensional processes where ]i5] ^ h? /{^h\,) the 
effective low-dimensional interaction can be analogously 
to the 3D case, replaced by a zero-range force. There- 
fore in each dimension D less than or equal to 3, a D- 
dimensional Bethe-Peierls approach can be relevant and 



useful as a tool for exploring the properties of few- and 
many-body shallow states. The dimensional reduction of 
the effective interaction from 3D to ID was first achieved 
in Ref. and from 3D to 2D in Ref. ^. 



II. ISOTROPIC CONTACT FORCES 

Two-body transition matrix and scattering 
amplitude 



One of the purpose of this paper is to show that the 
zero-range potential approach can be handled in a very 
simple and unified way in any dimension D < 3. The 
key tool used along these lines is introduced in Sec. [TT] 
of the paper. It consists of the general expression of the 
Bethe-Peierls asymptotic condition which can be used in 
any representation for all dimensions D. Substitution 
of the J-source terms of the zero-range approach by a 
family of well-behaved functions of vanishing but finite 
support allows one to avoid technical problems concern- 
ing double-limit calculations {i.e., evaluation of infinite 
series in the zero-range limit). While used usually in the 
configuration space, the zero-range approach can thus be 
implemented directly in the momentum representation 
for translation invariant systems. The link between this 
formulation of the Bethe-Peierls approach and the A po- 
tential introduced in Ref. f2q| is analyzed. The method 
is illustrated in Sec. IIIII with the example of the dimen- 
sional reduction issue from 3D toward ID and 2D. The 
Bethe-Peierls contact condition is directly expressed in 
the functional basis which diagonalizes the free Hamilto- 
nian and known results are obtained in a simple manner. 
For iV-body systems, integral equations which encapsu- 
late the Bethe-Peierls contact condition are obtained in 
Sec. IIVI It is also shown how to express this zero-range 
force approach in a second quantized representation. As 
an example, energy theorems [27h31| are derived by us- 
ing this last formulation. Section-fVl of the paper is de- 
voted to the three- and four-body shallow bound state 
issue. For three identical bosons in 3D, the original Skor- 
niakov Ter-Martirosian (STM) equation [S^l appears as a 
straightforward application of the results of the preceding 
section. A nodal condition is imposed on the three-body 
wave function in order to avoid the Thomas collapse [33| 
while keeping the simplicity of the zero-r ang e approach 
[3J]. Known results of universal theory [35| are recov- 
ered from the regularized STM equation. For two-mass 
component fermionic systems in 3D, critical mass ratios 
for the Efimov threshold are given in each partial wave 
by using the STM equation. Finally, the formalism is 
applied to derive four-body equations in the zero-range 
limit. In 2D, the eigenenergy equation for two-mass com- 
ponent tetramers is studied numerically. The energy of 
the ground s wave of '3-1-1' bosonic tetramers in the zero- 
range limit is computed as a function of the mass ratio. 



In this paper the norm of any vector v is denoted v. 
For a uni-dimensional system (D = 1), v is algebraic and 
ti = Ivl is the absolute value of v. 



This section reviews basic notions of two-body scat- 
tering in £)-dimensional spaces which are useful in the 
formulation of the Bethe-Peierls approach. 

Two colliding particles are described in their center- 
of-mass frame by the relative particle of reduced mass fi 
and relative momentum k. They are supposed to interact 
through a short range potential V, and the Hamiltonian 
for the relative particle is iJ = Ho + V, where Ho is the 
free Hamiltonian. In what follows, as a consequence of 
the small coUisional energy, scattering is supposed to only 
occur in the s wave channel of the relative particle. For 
an incoming atomic wave of momentum kg and of col- 
lisional energy E = h?kQ/2^ (measured from the contin- 
uum threshold), the scattering state |^ko) of the relative 
particle verifies the Lippmann-Schwinger equation: 



l*i< 



-|ko) 



1 



iO^ 



Ho 



-y|*i 



(2) 



{^\'^)-J-^i^\'^)^M^^-r). 



(3) 



We now turn to the k representation with the convention 
(r|k) = exp(ik • r): 

{2^Y 

In the k representation, the scattering states in Eq. ([2]) 
can be written as a function of the half on-shell transition 
matrix (or t-matrix) defined by: 

(k|T(ii; + jO+)|ko) = (k|F|vI/k„), (4) 



so that: 



(k|*ko) = (27r)^5(k-ko 



(k|T(£; + iO+)|ko) 



E + i{)+ - 



2fi 



(5) 



The pairwise interaction is supposed to be short 
range with a typical radius bo. For a small energy 
|i?| <C fi^ /{fJ,bjj) and momentum \li\bjj <C 1, the half on- 
shell i-matrix only depends on the collisional energy E 
and coincides with the low energy on-shell t-matrix de- 
noted as Tn: 



(k|T(S + iO+)\ko) ~ Td{E + tO^ 



(6) 



At the lowest order in energy (implying that |fco|6£) <C 1), 
the expression of Td can be parameterized by the D- 
dimensional scattering length ao- It can be written as 
|3l|33: 



TDiE + iO+) = 



0-3 



2/i 



X < 



(1 + ia3ko) 

-1 

lii{~ia2koeJ /2) 

-ikp 
^ (l + iaiko) 



{D = 3) 

P = 2) 



P = l), 



(7) 



where ^d is the full D-dimensional space angle: 
r^s = 47r, ^2 = 2tt, and fli = 2. 



(8) 



In Eq. ([7]) for the two-dimensional case, 7 is the Eu- 
ler's constant, and the 2D scattering length 02 is always 
positive. The transition matrices in Eq. ([7]) are the ba- 
sic objects for describing low energy scattering processes 
in D-dimensional few- and many-body systems. Scat- 
tering properties can be equivalently described with the 
_D-dimensional scattering amplitudes defined by: 



/2(fco) = l/lnHfcoa2eV2), 
/i(fco) = -l/(l + *fcoai). 



(9) 
(10) 

(11) 



The notion of scattering amplitude is often used in the 
configuration space where the scattering states for large 
relative coordinates (r » bo) are deduced from Eq. ([7]): 



(r|*i 



piko-r _ 



^iknr 



P = 3) 



foiko) X < 



r(i)r 



gifeor (23 = 1) 



(12) 

In Eq. (|12p , for £) = 2 the outgoing scattering wave func- 
tion Hq is the Hankel's function of order zero [3^; for 
D ^ 1 , fi is also called the even scattering amplitude [J] . 



B. Bethe-Peierls asymptotic approach in arbitrary 
representation 

Surprisingly the zero-range approximation has been 
essentially used in the configuration space while other 
representations can yield substantial simplifications. For 
example, as a consequence of translation invariance, the 
momentum representation is very well suited for solving 
the few-body problem in homogeneous space. In this sec- 
tion, it is shown how the zero-range approximation can 
be formulated in arbitrary representation. 

While for finite range forces Eq. ^ is only valid for 
a small relative momentum k {i.e., for kbjj <C 1), in the 
zero-range potential approximation this equality is ex- 
tended for arbitrary large values of k. That way, the 
short range pairwise potential V is replaced by a for- 
mal zero-range potential such that the scattering states 
in this approximation coincide with Eq. ^ in the small 
relative momentum limit (khr) <C 1) or with Eq. (J12p for 
large interparticle distances (r ^ bo). The simplest way 
to implement the zero-raiige approximation is to follow 
the Bethe-Peierls method [1|, where for any state j'l') de- 
scribed by the Hamiltonian Hq + V, the pairwise short 
range potential V is replaced by a delta-source term \d^) 
with an amplitude S^: 



V\^)^S^\d''), 



(13) 



and S'lj, is deduced from an asymptotic condition on the 
state I*). In Eq. ((13]), {S'^) is such that the usual S- 
distribution is obtained in the configuration space: 



{r\6^)^S{r), 



(14) 



and I J^) is also denoted by the formal ket of zero relative 
coordinates |0r). From Eq. ([5]) the scattering state |^ko) 
is written in the zero-range potential approach as 



l*ko) = |ko 



Td{E + iO+) 
E + ?:0+ - Ha 



15^). 



(15) 



Injecting Eq. 
tion ^ gives 



T3|) in the Lippmann-Schwinger equa- 



l*ko) = |ko) 



iSxJi 



ko 



E + M+ - Ho 



\S^) 



(16) 



Identifying Eq. (TTSl) and ((T6| shows that the source am- 
plitude associated with a scattering state |^ko) is nothing 
but the on-shell transition matrix: iSi^^ = Td{E + iO'^). 
In configuration space, scattering states in Eq. (|15p or 
in Eq. (1121) are singular for vanishingly small interpar- 
ticle distances. As shown in what follows, this singular 
behavior defines the Bethe-Peierls asymptotic condition. 
From Eq. (fT5)) one finds that all the scattering states for 
a given dimension D have the same singular behavior 
which reduces to 



(rl*) 



03 



Ax < 



In 



0-2 



-fO(r) {D = 3) 



-fO(r) {D = 2) 



(17) 



l(r-ai) + 0(r) (i^ - 1). 



In Eq. (fT7|) A is a function of the energy E but does not 
depend on r. It is important to note that Eq. (|17p is 
obtained in the center-of-mass frame, meaning that the 
limit r ^- is taken at fixed value of the center-of-mass 
of the two colliding particles. Conversely, for a given 
dimension D, Eq. ()17p is the sufficient asymptotic con- 
dition which permits one to obtain the source amplitude 
S^^ of a scattering state |^ko) at energy E in Eq. (fTBI) . 
The asymptotic condition in Eq. (flT)) referred also as the 
"contact condition" or the Bethe-Peierls condition, re- 
mains unchanged for any linear combination of scattering 
states and can be thus used to find any eigenstate in the 
zero-range approximation for a system where an external 
potential is included within the free Hamiltonian Hq. For 
a positive energy E, the Lippmann Schwinger equation 
is 



I*) = |*(")) 



0\]> 



E + iO+ - Hf 



■|'5^), 



(18) 



where |\['*-°)) is the complementary solution i.e., it is an 
eigenstate of the free Hamiltonian Hq at energy E (for a 



negative energy j^*^"^) = and the +iO~^ prescription in 
Eq. p^ is useless). The relation between the number A 
in Eq. (fT7| and the source amplitude 5* in Eq. ([T8|) is 
given by 



A = 






(19) 



The crucial point of this section is to show that the 
asymptotic condition in Eq. (IT71) can be expressed in ar- 
bitrary representation. In what follows for convenience, 
the delta distribution is represented by the e — >■ limit 
of a Gaussian weight: S^{r) — lime_>.o(r|(5f ) and 



(r|<5f 



(27re2)^/2 



exp 



(20) 



With this particular choice, in the momentum space the 
representation of the delta term has the same expression 
[denoted Xe{k)] for all dimensions; 



(k|<5f)=x.(fc)-cxp - 



(21) 



Matrix elements and states without index e are con- 
sidered in their zero-range limit which corresponds 
in this formalism to the limit e — > (for example 
lim(:_j.o |v&e) = I^P). For the formulation of the zero-range 
approximation it is also useful to introduce the 'reference 
state' denoted \(j)^) which results from the action of the 
two-body Green's function in free space on \dP) at the 
negative energy Ea: 



\^'. 



I 



Ea-Ho 



\6^) where Ea 



2/i 



<0. (22) 



In Eq. (|22|) the parameter A is chosen positive (A G n^ ; 
according to the usual prescription in scattering theory. 
In the momentum representation: 



(k|^^) - -^ 



2/i 



Xejk) 
fc2 -I- A2 ■ 



(23) 



In the zero-range limit where e — > 0, in the configuration 
space the reference state has the same type of singularity 
as inEq. pT]): 



2/i 



--+A]+0{r) {D = 3) 

r ' 



In 



Kre^ 



0{r) {D = 2) (24) 



r-j]+0{r) {D = l). 



The contact condition in Eq. (|17p can then be written in 
terms of the reference state as 

hm lim(r|*, - S*0^) = -^. (25) 



In the configuration space, the state |\l/e — S^icj)^) is a 
smooth function for (r, e) close to (0, 0). It is thus possi- 
ble to commute the r — > and e — >■ limits in Eq. (f25|) : 



lim lim(r|\I'c 



s*<^^) 



lim lim(r|^e 



S*</)^). (26) 



For particles moving in the presence of an external poten- 
tial, the source amplitude in configuration space depends 
on the center-of-mass coordinates of the interacting pair. 
Hence, in general situations the source amplitude is re- 
placed by a state associated with the center-of-mass of 
the interacting pair: 



D\J/ 



IS*) 



(27) 



Finally, the Bethe-Peierls condition can be written with- 
out specifying any representation as: 



lim(Or|*e 



S*(; 



IS*) 



ToiEA)' 



(28) 



Equation (|28|) is by construction, invariant in a change of 
A e M+ : as is shown in the next section, one recovers the 
so-called A freedom of the A-potential [2g|. This freedom 
permits to simplify exact calculations without introduc- 
ing any approximation, or also to improve approximate 
schemes [26|, |40|, |41|. Moreover, the e — > limit can be 
taken equivalent ly as follows. 



lim(5f|*, -S*<; 



IS*) 



Td{Ea)' 



(29) 



or also by substituting the ket in the left hand side of 
Eq. dill) by its e ^ limit: 



lim((5f|*-S*<; 



IS*) 



Td{Ea) 



(30) 



The contact conditions in Eqs. ([28l) - ([30)) can be expressed 
in any desired representation by inserting a closure rela- 
tion in the scalar product concerning the relative particle. 
For example in the momentum representation, one can 
insert the closure relation 



d^/t 



|k)(k| 



(27r)^ 
at the right of the bra {5^ \ in Eq. (PUI and this gives 



(31) 



lim 



{27t)D 



X.(fc)(k|*-S*<^^ 



IS*) 



Td{Ea 



(32) 



In the next sections, it is shown that Eqs. ([28)) - (|30| al- 
low one to obtain in a simple way the standard integral 
equations of few-body problems in the zero-range limit 
and also to compute the induced scattering resonances 
and related scattering problems in the presence of a har- 
monic transverse confinement. As a conclusion of this 
part, it is interesting (and useful) to note that for D = 3 
or D = 1, it is possible to eliminate the presence of the 



reference function in the contact condition by performing 
a specific limit on A. For I? = 3, in the A — >■ limit the 
contact condition can be written as: 



lim 9e (e((5, 

e-i-O ^ 



D=3 



*)) 



2Trh?a3 



(33) 



Equation (|33p can be written in the following alternative 
form: 






(34) 



where Reg extracts the regular part of ((5^~'^|^) in the 

e-J-O 

limit where e — >■ 0. For D = 1, elimination of the refer- 
ence function occurs in the limit where A tends to cxd. 
One obtains: 



lim(5f=i|vI/) = -^|5*). 



(35) 



C. Link with the s wave A potential 

A zero-range s wave pseudo potential was introduced 
by Fermi [42!] in order to perform calculations in the first- 
order Born approximation. Thus, the Fermi pseudo po- 
tential cannot be used for a non-perturbative approach 
(for large values of 03, for example). In the modern for- 
mulation of the zero-range pseudo potential, the Bethe- 
Peierls condition is included in the Schrodinger equation 
via the zero-range pseudo potential itself 0, Q ■ I* has 
been shown in the configuration space that there exists 
in each dimension D a family of zero-range potentials: 
the so-called A potentials .26,]. This section, links the 
Bethe-Peierls approach and the A potential in arbitrary 
representation. For this purpose, a family of operators 
R^ is introduced. They act on a state |^g) as: 



(36) 



i?^|vI/,) = |5f)lim((5f|vI.,-S*^ 

e-i-O 



For a state verifying the contact condition in Eq. ((29)) 
one has: 

i?(i/^°)|*,) =0 (37) 

where ao is related to the scattering length az) as 



ao 



ajj for Z? = 3 or D = 1 

a2e''/2 for £> = 2 



(38) 



These regularizing operator satisfy two other properties: 
(?) for a regular state (S** = 0) thus R^\^) = |(5f )(0r|*) 
and (m) i?^|0^) = 0. 

Explicit expressions for R^ can be found without spec- 
ifying any representation as follows: 



Rt\^.) = l-^f 



lim r, 

e-fO 



[(^f l*e)] , 



(39) 



where r^ is an operator defined in each dimension as: 



A £• 



= < 



^ ln(e^A2eV2)a, 



P = 3) 



{D^2) 



n 1 
2 A 



5. 



(D 



(40) 

Expressions of r^ in Eq. PO]) depend on the choice 
made for the short range functions (rjjf'), given here by 
Eq. (pO)) . The source term in the stationary Schrodinger 
equation can be expressed in terms of R^ by using 
Eqs. dm and [Ml): 



{Ho - E) I*,) + Tz,(^A)i?f |*e) = 0. 



(41) 



In the limit e — >■ 0, the pseudo potential in this equa- 
tion coincides exactly with the A-potential introduced in 
Ref. [26[, where the coupling constant (^a is nothing but 
the transition matrix at energy £'a: 



3A =Td(Ejs 



D. Context of resonant scattering 



(42) 



The zero-range approximation is especially interesting 
in the regime of parameters where the pairwise potential 
leads to a resonant scattering at low collisional energy. In 
this regime, the cross section proportional to |/dP has 
a large value at small relative momentum {kobjj <^ 1) 
and expressions in Eqs. ([9])- (fT2)) are good approxima- 
tions. For Z? = 3, this corresponds to the regime where 
the scattering length as is large in absolute value as com- 
pared to the potential radius 63 and the maximum is at 
ko = 0. This regime can be achieved by using the FR 
mechanism (see Ref. [20] for a general review of FR in 
ultracold atoms) . The FR involves the coupling between 
atoms (in the "open channel" ) and a molecular state (in 
the "closed channel" ) characterized by a size of the order 
of the potential range 63. This two-channel description 
provides the expression of the scattering length (03) as 
a function of the external magnetic field B in the vicin- 
ity of a given resonance located a.t B — Bq which can be 
parameterized by the formula [43| 



03 = abg 1 



_AB_ 
B-Br 



(43) 



In Eq. (|43p , AB is the width of the magnetic FR and Obg 
is the "background" scattering length i.e., the scattering 
length away from the FR. The magnetic width (AB) can 
be also characterized by the "width radius" R* defined 
by 



R* ^ 



ma]jgSfj,AB ' 



(44) 



where S^jl is the difference of magnetic moment for an 
atomic pair in the open- and the closed-channel [20, El] . 
In the vicinity of a broad resonance, R* is of the order of 
(or smaller than) 63 the 3D scattering amplitude can be 
approximated by Eq. (|9l) in a large interval of momentum 
where |fco|53 ^ 1. 

As shown in the next section, the resonant behavior 
in the low dimensional atomic waveguide can be reached 
by tuning the length fltrans of the transverse confinement 
and/or the 3D scattering length. For D = 1 the reso- 
nant regime occurs in the limit where the ID-scattering 
length ai is small with respect to the transverse length 
(i.e., I Oil ^ 61) and also for a vanishing relative momen- 
tum fco — > 0. In 2D, the maximum in I/2P occurs at 
the momentum fco = 2/(e^a2) (where I/2I = 2/7r) and the 
resonant regime is thus reached at low energy if and only 
if the 2D scattering length is large as compared to the 
transverse length (i.e., 02 3> b2). 

In the resonant regimes defined above, the probability 
of scattering for two colliding particles is large so that 
the pairwise short range potential affects the form of the 
wave function at interparticle distances which are large 
with respect to the potential radius bi). More precisely, 
the scattering states in Eq. ([5]) can be decomposed in an 
incoming part |0inc) = |ko) and a scattered part |0scatt)- 
At large distance and in the resonant regime: 



(r|0i, 



r»fc£ 



laal/r {D ^ 2,) 

Vi^/V^ {D^2) (45) 

1 {D^l). 



For D = 3 or D = 2 this ratio is greater than or of the or- 
der of unity for feu ^ J' < | a_D | , where | od | is arbitrarily 
large and is equal to unity for arbitrary large distance in 
one-dimensional systems {D = 1). The zero-range poten- 
tial approach is a formalism which permits one to evalu- 
ate accurately the wave function in configurations where 
particles are outside the potential radius while configu- 
ration where two or more particles are inside the poten- 
tial radius are not reliably described. Consequently, the 
Bethe-Peierls asymptotic approach is very well suited for 
studying systems in the resonant regime where the wave- 
function is modified by the interaction at interparticle 
distances which are large as compared to the potential 
radius. 

In the 3D space and for a positive scattering length, 
the on-shell transition matrix T3 in Eq. ([T]) has a real 
pole at negative energy E2 — ~h? /{^^a^). In the reso- 
nant regime, this pole is associated with the existence 
of a shallow dimer which is thus very well described in 
the zero-range approach (the probability that the rel- 
ative pair has a radius greater than the potential ra- 
dius &3 is exp(— 263/03) ^ 1). For D — 2, the pole is at 
E2 ~ — 2?i^/(/ia2 exp(27)), and in the resonant regime it 
is also the signature of the existence of a shallow dimer. 
For _D = 1, in the resonant regime (oi — > 0) the approxi- 
mation of the transition matrix in Eq. ([7]) doesn't possess 
a negative low-energy pole and next-order terms in the 



energy expansion are needed to find the pole correspond- 
ing to the lowest bound state [45{ . 



III. DIMENSIONAL REDUCTION 

In this section, the dimensional reduction issue is 
solved by using the zero-range approximation with the 
Bethe-Peierls condition expressed through Eq. ([M)) . This 
section gives alternative derivations of the results given 
in Refs. 0-01 • The problematic of the 3D^1D (respec- 
tively 3D— )-2D) dimensional reduction is as follows: two 
particles move freely in ID (respectively 2D) while they 
are confined in the transverse direction by an atomic 
waveguide built from a 2D (respectively ID) trapping 
potential. The collisional energy is such that at large 
interatomic distances, the particles are confined in the 
ground state of the trapping potential. Solving the issue 
amounts to finding the low energy transition matrix of 
this low dimensional quasi-2D (resp quasi-lD) scattering 
process as a function of the waveguide parameters and 
of the 3D scattering length 03. In what follows only the 
case of harmonic trapping where the center-of-mass and 
relative motions decouple is considered. In the center-of- 
mass frame the source amplitude ^^r is constant and the 
scattering state can be written as: 



I*) = |*(o)) + 1*"^ 



(46) 



where j^E''*'^) is the incoming wave and the "interacting 
part" is 



I* 



int\ 



Si^ 



E + iQ+ -Ho' 



lOr 



(47) 



In Eq. (j47p . the source amplitude S* is found from the 
Bethe-Peierls condition by solving the equation: 



(Or|*(°') 



Regis, 

e-i-O 



D=3 



* 



int\ 



271^2 



as 



(48) 



These quasi-lD and quasi-2D scattering problems have 
been initially solved by using the Bethe-Peierls method 
in the configuration space in Ref. [4] and Refs. [a, 0- In 
this section, an alternative derivation is given in order to 
illustrate the formalism introduced in the previous sec- 
tion. 



A. Linear atomic waveguide 

The 3D— >-lD reduction problem, is considered here in 
the case where the two colliding atoms are confined in 
an isotropic two-dimensional harmonic trap while they 
move freely along the third direction (z). In this case, 
the problem can be solved in the center-of-mass frame 
where the noninteracting Hamiltonian is 



Ho 



2jLt 



52+H2D, 



(49) 



and H2D is the Hamiltonian for the transverse motion: 
H2D = -^{dl + dl) + \^lLol{x^ + y^) - hLo^. (50) 

The hnear atomic waveguide is characterized by the 
transverse length: 



a_L 



/ia;_L 



(51) 



and the zero-range approximation of the D ~ i pair- 
wise potential is justified in the limit where a_\^ 3> ^3. 
In Eq. ([50)1 . the zero-point energy has been subtracted. 
That way, the energy E of the scattering states is mea- 
sured with respect to the continuum threshold. In this 
quasi-lD scattering problem at energy E = h^k^/2iJL, the 
incoming wave for the relative particle (state j^'"^)) in 
Eq. (|46p has a momentum kg along z and is in the ground- 
state of the transverse 2D oscillator: 



|*(0)) = |n, = 0)|% = 0)|fc, = fco), 



(52) 



where Ux and \nx) (respectively Uy and l^j^)) are the 
quantum number and eigenstate of the ID harmonic os- 
cillator of atomic frequency a;_L and mass /i along x (re- 
spectively y). The system is quasi-lD for collisional en- 
ergies such that: 



E<E, 



ID 



(53) 



where El^^^ — 2haj±. Equation (|55|) defines the 
monomode regime of the atomic waveguide such that out- 
going particles are trapped at large relative distances in 
the ground state of the transverse trap. The contribution 
of the incoming wave in Eq. (|48p is: 



(Orl* 



(0)\ _ 



0o(O)2 = 



1 



y/7m± 



where ipnix) — {x\n) is given by 
exp 



0™ {x) 






(54) 



(55) 



In the calculation of (jp^^j^I''"*), the quantum numbers 
of the noninteracting Hamiltonian are introduced directly 
by inserting the closure relation: 



+00 



27r 



%){K\ J2 l"-)KI E l"«)KI=I- (56) 



Only even values of Ux and Uy give a non vanishing con- 
tribution, and one obtains 



(S, 



D— 3 l^int 



hujj_ J _ 



+00 



dkz sr^ 

p,q=0 
cD=l 



exp 



fc?e2 



(5f=i|2p)(J,"=i|2q)</>^p(O)0^,(O) 



fc2Q2 



(57) 



where the dimensionless energy variable r is defined by 

E 



2huj± 



iO^ 



(58) 



The behavior of (jp^^j^I''"') in the limit where e van- 
ishes is more easily obtained in the domain of negative 
collisional energies (r < 0) than in the domain of positive 
collisional energies E > and are related the each to the 
other by analyticity. For r < one can use the identity 



fc5o2 



p 



due 



^+n+p—T 



(59) 



which allows one to decouple the discrete summations 
from the integration over k^ in Eq. (j57p . From Eq. (j55l) , 
one can deduce the following limit: 



limV(5?=i|2p)0;^(O)e 



-pu 



\M0)\' 



(60) 



Using Eq. (IM)) and integrating over kz in Eq. ([57]) gives: 
— /iSij, f°° du 



{S. 



D— 3 lij/int 



^0 27r3/2ft2 



a± Jo 



1 



(61) 

In the limit where e — > 0, this expression diverges as 1/e 
and its regular part coincides with a Hadamard's finite 
part 46, 47]: 

mS* „ „ f°° du 



Regis, 

e-i-O 



D=3| 



xpintN ^ 



27r^/^h^ai_ 



p.f. 



/ul 



(62) 

One recognizes in Eq. (|62p an expression similar to the in- 
tegral representation of the Hurwitz zeta function ^(s, z) 
defined for 3?(s) > 1 and < z < 1 by H: 



C(5,^) 



r(.) 



dt- 



1-e- 



(63) 



However, in Eq. (|551) the integral diverges for s = 1/2. 
The Hadamard's finite part in Eq. ([S^ permits one to 
achieve a meromorphic continuation in the variable s 
of Eq. §^ and to identify the function C(l/2,^) [ioj . 
Hence, 



S\]> — 



2^i/2;j2 



c(i-T) + ft 



(64) 



For positive energies and in the monomode regime 
{0 < E < 2haj±), the scattering state written in the con- 
figuration representation in Eq. (|46p has a non evanescent 
contribution at large relative length z ^ a±, in the sub- 
space of the transverse ground-state only. This property 
allows one to identify a quasi-lD scattering process with: 



{z^Hx = 0,ny 



0|*) = exp(ifco2;)- 

\z\:»a± 



z/iSij. 



y/T^fi^ax^kQ 



exp(ifco|z|), (65) 
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and a± plays the role of a ID-potential radius {i.e., 
bi ^ a_i_). The low-energy scattering amplitude j3D-j-id 
of this quasi-lD system is thus: 



/ 



3D->1D 



(ko) 



z/iS.1. 



^y^^h?a±ko 
In the |A;2)-representation, the scattering states read 

{k,,n^ = 0,ny^ 0|*) ^ {2n)S{k, - fco)+ 

Sv]/ 



(66) 



y^a± (e + iO+ 



2m 



(67) 



so that the quasi-lD scattering transition matrix y3D->-iD 
can be expressed as 



nSD-i-lD 



S\]> 



/7ra_L 



Finally, from Eqs. ([M)) and (|68p one obtains: 



T 



3D^1D 



2^2 



A^«± C(i-r) + f^ 



(68) 
(69) 

(70) 

and for positive energies ^^t = —ikQa±/2. From 
Eq. ([7]), one can then identify a ID-scattering length 
which is a function of the transverse length a± and of 
the 3D scattering length 03 ^4^ : 



In the low energy limit {koa± <C 1), 



^[I'-^Jr^oW-r 



'-a 



ai = 



a± 



as V2 



(71) 



B. Planar atomic waveguide 



\nz) is the eigenstate of quantum number riz for the ID 
harmonic oscillator of frequency w^ and mass fi. The 
subsequent derivation supposes that the 3D zero-range 
approximation is justified, i.e., that the atomic waveg- 
uide is such that: 



dz > ^3, 



(74) 




where Gz is the characteristic length of the atomic waveg- 
uide: 



(75) 

In the quasi-2D scattering problem the incoming state is: 

|*(°')-|k2D-ko)|n, = 0), (76) 

where kg is the 2D- wave relative wavevector of the incom- 
ing wave and k2D = k^e^ + kySy is the 2D collisional rel- 
ative momentum. The collisional energy is i<^ = fi^k'^/2^ 



and the monomode regime condition is given by: 



E <K 



2D 



(77) 



where i^™ns = ^z. Equation (|77l) ensures that the sys- 
tem is quasi-2D; i.e., the colliding particles are asymp- 
totically trapped in the ground state of the transverse 
trap. The behavior of ((5P"^|^'"*) in the limit where 
e — ^ is obtained by using the same techniques as for 
the 3D— >1D reduction problem. The quantum numbers 
of the free Hamiltonian are introduced by insertion of the 
closure relation: 



«2D 



(27r) 



k2D)(k2D| 



\nz){nz\ 



(78) 



In the 3D Bethe-Peierls asymptotic condition, only even 
values of riz contributes and one obtains: 



In the 3Dh>2D reduction problem, the two colliding 
particles are confined in a planar harmonic waveguide 
of frequency Uz along the z direction while they move 
freely in the two other directions. The noninteracting 
Hamiltonian in the center-of-mass frame reads: 



Ha 



2n 



id'. 



dl) 



■Hz, 



(72) 



with, the ID-harmonic trap along z. 



Hz 



2^l ' 2 "" 



'2 



-^- (73) 



In Eq. ([75]) the zero point energy of the transverse trap 
has been subtracted, so that similarly to the previous 
3D— >1D reduction problem, the energy of an eigenstate 
is measured with respect to the continuum threshold and 
coincides with the collisional energy. In what follows. 



(S, 



D— 3|ayint 






X ^exp 

p=0 



^^20^ 



fc2D&2D 

(<5P=i|2p)</)! (0) 



(79) 



P 



where the dimensionless energy variable r is defined by 

E 



2hUJy 



iO^ 



(80) 



The discrete summation in Eq. (j79p is performed in the 
domain r < by using the transformation of Eq. ((59)) 
together with the identity of Eq. ([60]). After integration 
over fc2D, the regular part of ((5g|^'"*) can be expressed 
in the limit where e — >■ as 



Reg(<5,3|*'"')=-|^|<Ao(0)|V(r), 



(81) 



where the function J{t) is defined in the domain t < 
by: 



J(r)=P.f.r 
Jo 



^du exp(r'u) 



u yj\ - exp(-M) 
From Eq. (|48p . one finally obtains: 



^ + |</.o(0)pj(r) 



(82) 



(83) 



Projection of the wavefunction on the ground-state of the 
ID transverse harmonic oscillator (k2D,"-z = 0|^) gives 
the quasi-2D transition matrix: 



T 



3D-i.2D 



7r?i2 



^ 



2a3 



J(t 



(84) 



In the domain of positive energy Eq. (j80p is continued 
analytically from the identity 



J{r) = In - 



B 

"hrr 



El- 



n \{2n- 1)!! 
(2n)!! ' 



(85) 



where B ~ 0.9049 is defined by f3£ 



In 



Bel 

'2^ 



du 



Vi 



i3/2 



1 



(86) 

and 7 = 0.5772 ... is the Euler's constant. The 2D low 
energy condition is given by: 



\E\ < hhJ:, 



(87) 



and az plays the role of a 2D-potential radius, i.e., 
&2 ^ Oz- Only the first logarithmic term in Eq. ([85]) con- 
tributes at the lowest order in the limit of Eq. (|87l) , and 
one can identify the resulting expression of the quasi-2D 
transition matrix X'3D->-2D .^j^^Jj^ ^^jj^g 2D low energy transi- 
tion matrix T2 in Eq. ([7]). The 2D scattering length 02 of 
this quasi-2D scattering problem can thus be expressed 
in terms of a^ and of the 3D scattering length 03 [y, [3, |41| 
as: 



02 




2a3 



(88) 



It is interesting to discuss a little bit further the con- 
dition for having a low energy resonant behavior in this 
quasi-2D system. The 2D cross- 'section' (this is a length) 
is proportional to |r'^^^^^(fco)P/fco and a maximum oc- 
curs in the low-energy regime Eq. (|87p for: 



a 



03 y/ii \BhWz 



(89) 



The right hand side of Eq. (|89| is large and negative 
which shows that the resonance occurs for sufficiently 
large values of |a3| with respect to a^ in the domain of 



negative 3D scattering lengths (03 < 0) . At fixed value of 
the coUisional energy E and of the 3D scattering length, 
the resonance can be reached by tuning the harmonic 
frequency of the atomic waveguide: this is the so-called 
quasi-2D confinement induced resonance first found in 
Ref. [|. 



IV. FORMALISM FOR iV-BODY SYSTEMS 

A. k representation 

In this section the Bethe-Peierls method is applied to 
few- and many-body systems in the k representation. 
The system is composed of N particles of respective 
masses mi,m2, ■ ■ .misr and of momenta ki,k2, . . . kjv. 
For convenience the set of momenta is denoted by the 
short-hand notation {k}. The relative and total mo- 
menta of a given pair of particles (ij), are denoted by 






-!»■ ij — Ki 



^J' 



(90) 



and the reduced mass of the pair (ij) is denoted ^ij 

mirrij 



(91) 



It is also convenient to introduce the set of momenta 
denoted {C} with {i'^} = CI ,Ci 



>2 



■ • Cn where: 



K, 



K, 



2 (mi 



2 (mi 



for n = i 

for n = j 
otherwise 



(92) 



This last notation permits one to isolate in the kinetic 
energy the contribution of the relative particle (ij) from 
the A^ — 1 other momenta: 



N 



n=l 






N 

E 

n=l 



m' 



(93) 



A non vanishing source amplitude |Sj[~'"') is associated 
with any pair of particles interacting via the short range 
pairwise potential in the many-body state |^e). In the 
case where the system is composed of particles with spin, 
in what follows |^c) denotes the projection over a given 
spin configuration of the many-body state including the 
symmetry imposed by the quantum statistics. The free 
Hamiltonian of the system (which may include an exter- 
nal potential) is denoted Hq, and the A^-body stationary 
state at energy E verifies: 



iHo-E)\^,)=~Y.\^ij) 



D\|Ql^, 



s^l 



sr)- 



(94) 



i<j 



In Eq. dHD, \{ij) : <5f ) is the state |(5f ) for the relative 
particle (ij) of momentum kij, and the source amplitude 
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|S^ -') is a state for the A^ — 1 other particles of momenta 

{C''}- If the pair (ij) does not interact then jS^-') = 0, 
otherwise it satisfies the Bethe-Peierls contact condition: 



hm 



i((zj):'5f|(|*)-|fe-):0^)|Sr^")) 






(95) 

where the reference state \{ij) : (j)^) is defined in the k 
representation as: 



(fc,,|(zj):0^) = - 



1^ ^ 






A2 



(96) 



In what follows Eqs. ([M)) and (1^51) are used in the case 
where there is no external potential and l^*^"^) denotes 
the complementary solution of Eq. (|94| which satisfies 
the boundary conditions of the problem considered (for 
example an incoming plane-wave for a scattering prob- 
lem). The many-body wavefunction can then be written 



(ij) dS^-')). By analytical continuation, the same inte- 
gral equation holds for i? > (in which case |\l/(°^) ^ 0). 
Finally, assuming that the e — > limit is well defined, for 
each interacting pair one obtains the following integral 
equation: 

E J ^SE{{n){{e'}\sm 



n<p 



yT|c«^J\ 



({^lis^* 



TDiE^^+^0+) 



{2^Y 



.({k}|*(°)). (101) 



One has to notice that the integrals in the first line of 
Eq. (|10ip are performed with the constraint that j^*-'} is 
held fixed (but {C^} is not !). 



B. Second quantization 



({k}|v[/,) = ({k}|vi/(o)) + c;^({k})Ex.(%)({r^}|s**^ 



-3\ 



i<j 



(97) 

where tj^({k}) is the A^-body Green's function in the k 
representation: 



rj + lU Z^n=l 2m„ 



QEim) = 



For each interacting pair, the contact condition Eq. ((95)) 
gives in the k representation an integral equation in terms 
of the source amplitudes and of the complementary solu- 
tion 1^^°^). Interestingly, one can extract the contribu- 
tion of the interacting pair {ij) in the contact condition 
Eq. (j95p without any integration by using the fact that 
this equation is invariant in a change of A G K"*" (i.e., by 
using the A freedom). To this end, one expresses the 
Green's function in Eq. (P7| as: 



Qe{{\^}) = 



1 



^^^^-~Wi 



(99) 



,fe) 



where E),^^ is the collisional energy of the pair (ij) de- 
fined by: 



E^=E-Y. 



"" f^m 



2mr, 



(100) 



It is the energy of the pair (ij) in its own center-of-mass 
frame while the other particles do not interact for a given 
set of momenta {k} and a total energy E. 

For a negative energy {E < 0) in which case l^*^")) = 0, 
without loss of generality one can make the particular 
choice A2 = -2^i,^.jE''^^l/h^ in Eq. ^. That way, the 
term involving the reference state (p exactly cancels with 
the term associated with the source amplitude of the pair 



In this section only fermionic particles of the same mass 
m and two spin-components are considered. General- 
ization to other systems gives similar equations. The 
creation and annihilation operator of an atomic wave 
of momentum k for a spin a are denoted a\^ ^ and ak.o- 
with the standard anti-commutation rule: {ak,o-, <^k ^'1 ~ 
(27r)^(5(k - k')(5crcr' and {ak,^, ak.o-'} = 0. 

The general expression of the A-potential is a simple 
way to write the Hamiltonian in its second quantized 
form. Using the definition of the A potential for finite 
values of e in Eq. (PT|) , the Hamiltonian can be expressed 
as: 



H,^ 



d"k 



E 



ti'e 



«k.<T«k,<T 



Td{Ea) f d'^K 



X lim r^ 



(27r)'f^ ^ 2m 

(27r)^^'^ JaiK-k',t"iK+k',; 

7^— TDXelKjOiK-k^i^iK+k.t ' 



(27r) 



(102) 



where the dot (•) after the annihilation operator a i j^_|_j(. ^ 
reminds that the limit e — ^ depends on the many-body 
state on which H^ applies. 

In the hypothesis where the few- or many-body state 
l^e) is well defined in the zero-range limit (where only the 
Bethe-Peierls contact condition is used), the mean energy 
of the system (\['e|7JJ_^f) is a regular function of e and 
the energy theorem |27H3ll | follows from the identity: 



lim r 



(1/5b) 



(vl/Ji/JvI/ \1 = IH) 



(103) 



The state |^e) verifies the contact condition for each in- 
teracting pair thus using the property in Eq. p7|) , action 



of the operator lime_j.o r, 



(i/ao) 



on the interacting term 
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in (^e|7Je|^c) gives exactly zero and: 



(H) = linir(i/"°) 



d^k .^ h^P I 



{2tt)D 2^ 2m ^"k..ak..;e 



(104) 

The energy of the system can be thus expressed in terms 
of the one-body density: 



nk,cr = hni(a! ak,cr 



(105) 



In the hmit where e — >■ 0, one obtains from Eq. (1^71) the 
high-momentum behavior of nk.o-: 



nv. 



fc— >oo h, 



o 



(106) 



where S is spin independent since the interaction only 
occurs between particles of different spin. Hence, the 
integrand in Eq. (|104p behaves as p- at high momentum 
and the action of the regularizing operator gives: 



{H) 



d°k 
(27r)^ 



E 



2m 



rik., 



h^Sdl 



2m{l + k^al,) 



(107) 
where the use of the length ao defined in Eq. (pS]) permits 
one to gather in a single form the energy theorems for the 
three different dimensions. 



V. FEW-BODY BOUND STATES 

A. Efimov states 

In this section, the generic case where the Efimov ef- 
fect occurs is considered: three identical bosons of mass 
m interact in 3D with a pairwise interaction of vanishing 
range with respect to the 3D scattering length 03 |50l.l5ll). 
These states have been observed for the first time in 
ultra-cold atoms |52h55|. In what follows, some univer- 
sal properties of Refs. [3^, [50, [5l| are recovered from the 
STM e quat ion by using the nodal condition introduced 
in Ref. [3J]. However, it is worth pointing out that de- 
viations from universal predictions which are observed 
in experiments can be taken into account through finite 
range models |56l - [59| . The energy of a trimer E — E3 
in its center-of-mass frame is negative (and less than the 
dimer's energy E2 if it exists) so that the complementary 
solution in Eq. pUTj) is zero (|^(°)) = 0). The binding 
wavenumber g > of a trimer is defined from: 



2„2 



E. 



K'q 



(108) 



The Bose statistics imposes that the source amplitudes 
for each pair of particle coincide with the same function: 



aSi |q1±^3 



23i |q2^3\ 



{{e'}K^') = {{e'wri = ({rnisr^). (io9) 



As a consequence of translation invariance, the source 
amplitudes in the center-of-mass frame can be written 
as: 

({€''}|Sr') - (2^)^<5(ki + k2 + k3) F(k3). (110) 
The integral eigenequation Eq. (|101|) takes the form: 



F(k) 



/3(n/'7^ + ^ 



= 87r 



F(u) 



dh 



(27r)3 u2 + fc2+k.u + g2- 



(111) 

Equation (fTTTj) is the so-called STM equation [32]. This 
equation is rotationally invariant and can thus be studied 
in each momentum sector. In Ref. 60], Danilov showed 
that as it stands Eq. (|llip is ill defined: it supports a con- 
tinuum of negative energy solutions in the s wave sector 
of F{k). Hence, the Bethe-Peierls asymptotic method 
which is at first sight adapted for modeling the three- 
boson resonant problem does not permit one to derive 
a well-defined eigenequation. The s wave component of 
the source amplitude F(k) is denoted as. 



/S^« 



0(fc), 



and 4>{k) verifies the integral equation. 



0(fc) 



h 



^2+3|i 



— / du(j){u)JC„(k,u), 
^ Jo 



(112) 



(113) 



where the s wave kernel of the STM equation is given by: 

' u^ + fc2 + q"^ + ku^ 



JCq{k,u) = yln 



''^ ^ ' k \ u'^ + k'^ + q'^ — ku 



(114) 



The fact that there exists a continuum of negative en- 
ergy solutions means that the Bethe-Peierls model is not 
self-adjoint for the three-boson problem. Danilov found 
a way to restore the self-adjointness by introducing a 
supplementary condition on the high momentum asymp- 
totic behavior of the function 4>{k) [60l- The method of 
Danilov is based on the fact that for all values of q and 
a 7^ 0, Eq. (|113p supports solutions with the asymptotic 
behavior: 



0(fc) 



Ak 



isQ — 2 



Bk 



-isQ—2 



(115) 



where A and B are two constants and sq solves the equa- 
tion sin(7rso/6) = VSsq cos(7rso/2)/8 {sq ^ 1.00624. . .). 
The zero-range approach is made self-adjoint if one fixes 
the asymptotic phase shift between the two conjugate be- 
havior k^''^" for all values of a and q. However, Minlos 
and Fadeev showed that even with this supplementary 
phase-shift condition the spectrum is not bounded from 
below [61| : this is the so-called Thomas collapse which is 
characteristic of zero-range forces [33|. In 1970, Efimov 
solved the three-boson problem in the resonant regime by 
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FIG. 1: Shape of the expession of k'^(t){k) in Eq. (|118p plotted 
in a semi-log plot as a function of k/q. 



introducing the notion of a three-body parameter which 
is also in the zero-range limit, a way to fix the asymp- 
totic phase-shift in Eq. ()115p . In the k representation it 
is referred in what follows as n* and is defined through 
the asymptotic behavior of (/)(fc): 



(j){k) ex —J sin 

fe— s-oo k 



So In 



k^/?r 



(116) 



By construction the three-body parameter is not unique 
[Eq. (|116l) is invariant in a change k* — >■ k* exp(7r/so)]. 
At unitarity (i.e., |a| — >■ cxi and 63 —^ 0), Efimov showed 
that the spectrum of trimers is characterized by an accu- 
mulation point at zero energy: binding wavenumbers of 
the trimers are related the each to the others by a scaling 
factor. 



9pe 



-(n-p)TT/so 



(117) 



where (n,p) e 1? . Interestingly, in Ref. [621 the analyti- 
cal expression for the source amplitude of the trimers has 
been found at unitarity: 



<P{k)^ 



9^ + ^ 



: sm 



So arcsinh 



29 , 



(118) 



In particular, this result shows that the choice made 
for the definition of k* in Eq. (|116l) is such that 
the spectrum of the zero-range theory at unitarity is 
qn = K* exp(n7r/so) with n dZ. The shape of the func- 
tion k'^(j){k) of Eq. (|118p is given in Fig. ([T]). For increas- 
ing values of k starting from fc = 0, the first zero of (t){k) 
is located at the momentum k — — ^ sinh(7r/so). Using, 
this property, it is possible to filter a set of solutions of 
Eq. (irn)) satisfying Eqs. PH)) and (ITT7)) for sufficiently 
high quantum numbers (shallow states) and with a spec- 
trum bounded from below. This is done by imposing a 
nodal condition on the eigenstates. 



0(fcr) = 0, 



(119) 



where the node ¥^^ is chosen among the set: 



kl^ 






peZ. 



(120) 



The position of the node in Eq. (I120p fixes the min- 
imum energy of the spectrum. For example if one 
chooses the nodal condition for p — 2, using the fact 
that exp(7r/so) — 2sinh(7r/so), the minimum energy Eq 
is almost equal to ~h?K* ^ exp(27r/so)/m [with a relative 
error 2exp(— 27r/so) — 3.8 x 10^'^) and the spectrum at 
unitarity is asymptotically (i.e., for large values of n] 
given by 



En = exp 



-2(?i - l)7r 



So 



,ne N. 



(121) 



For a finite scattering length 03 and for fc|a3| ^ 1, the 
eigenfunctions 4>{k) of Eq. (|113p have the same behavior 
as the unitary solutions of Eq. (|118|) so that Eq. (I119P can 
be also used as a filtering condition and permits one to 
recover the universal spectrum of the zero-range theory 
for energies much larger than Eq. In principle, the zero- 
range theory also called 'universal theory' in the litera- 
ture is recovered by imposing the filtering condition at an 
arbitrary large node k'^J'^. To summarize this discussion, 
the nodal condition in Eq. (J119I) has two roles: first, it 
imposes the Danilov's asymptotic phase shift for states of 
sufficiently high quantum number and second, it imposes 
a minimum energy to the spectrum. For a realistic finite 
range force, the possible values of the binding wavenum- 
bers obtained from the zero-range theory in Eq. (|117p are 
such that 5^63 <C 1, the ground state energy in Eq. (|12ip 
has thus to be chosen higher than —h^/{mh\). 

The nodal condition in Eq. (I119P applied in Eq. (|113p 
can be transformed into an integral condition: 



2 r°° 

= / du(l){u)lCq{kl^«',u). 

T^ Jo 



(122) 



Subtracting Eg. (| 122l) from Eq. (I113P gives a regularized 
STM equation IjP 

0(fc) 



3fcf 
4 



2 f°° 

- / du [lCg{k, u) - lCq{kl°^, u)] (t){u), (123) 
■"■ Jo 

which encapsulates the nodal condition. The Danilov- 
Efimov contact condition (I116P is implemented exactly 
in the STM equation in the limit where the integer p 
tends to infinity in Eq. ()123p . Numerical solutions of 
Eq. (|123p are obtained by introducing an ultraviolet cut- 
off Q in the integral. One can verify that results are in- 
sensitive to the choice made on the cut-off for Q ^ kp^^. 
For instance considering the unitary limit, taking the 
nodal condition at A:™^ and an uv cut-off at 5 x 10^ x k*, 
one finds the relative error for the ground state equals 
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to 3.8 X 10"'^ ~ 2exp(— 27r/so) and for the first excited 
states (n < 8) one finds a relative error less than 10""'. In 
practice, the value of the scaling factor exp(27r/so) ~ 515 
is relatively large so that for the nodal condition at A:2°^i 
the spectrum is very close to the universal spectrum be- 
ginning from the second branch (n > 1) and for an inverse 
scattering length l/|a3| smaller or of the order of k* . 

In Fig. ([2]), the second Efimov branch of the regular- 
ized STM Eq. (fT23l) with the nodal condition at ^ has 
been plotted as a function of (l/as). As in Ref. [33|, the 
thresholds for the appearance of the trimers are denoted 
by a* at the atom-dimer continuum limit and by a' * 
at the three-atom continuum limit. As a consequence 
of the choice of the nodal condition at fc™s for a finite 
p, results slightly differ from the zero-range theory in 
Ref. ^35|. The atom-dimer scattering problem (regime 
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1/a^ ( units of k* ) 



Using the regularized STM equation one obtains 



3g{k)a3 



1- 

°° du 

■KV? 



{ka^y 



1 + (fca3)2 1 + {kl'^^a^y 



[/Ci/,3 {k, u) - /Ci/,3 {k;^^ u)] g{u). (125) 



The atom-dimer scattering length denoted (oad) is the 
zero momentum limit of the function g{k) and is plot- 
ted in Fig. ^ for the nodal condition taken at k'^^. 
Figure ([3]) is limited to atomic scattering lengths as 
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FIG. 3: Atom-dimer scattering length Had computed with the 
nodal condition at k^^ and plotted in a semi-log scale as a 
function of k* ja-j, for a^k^^ <^ 1. 



FIG. 2; Solid line: Second branch of the trimer's spectrum 
obtained with the nodal condition at ^2"^. Vertical dashed 
line: Trimer's appearance thresholds. Oblique dashed line: 
Atom-dimer continuum limit. The value of the scattering 
length at the trimer's appearance threshold a* ~ — 1.51/k* 
and a' * ~ .0713/k* are close to the results of the "universal 
theory" in Refs. ^,\^, a* ~ -1.507/k* and a' * ~ .0707/k*. 
The difference is due to the fact that the zero-range theory 
is recovered only for a nodal condition at an arbitrarily large 

l.i-eg 



where 03 > 0) can be also solved by using the same fil- 
tering technique. The atom-dimer scattering length is 
computed as a function of the atomic scattering length 
03 for a given value of the three-body parameter k* . To 
this end, the atom-dimer collisional energy is set to zero 
{q = 1/03) and the ansatz for the source amplitude is [32l| 



0(fc) = 27r 



2'5(k) 



fc2 



-An 



9{k) 



(124) 



much smaller than l/k2 i.e., in a regime where the 
nodal condition permits to recover the universal the- 
ory with high accuracy. The scattering length Oad di- 
verges at the threshold of appearance of an Efimov's 
trimer [a3 ~ a* exp(n7r/so), n € N] and thus exhibits the 
log-periodicity which is a characteristics of the Efimov 
physics. 

In Ref. [63] another regularizing technique of the STM 
equation was derived in the framework of the effective 
field theory. In this last formulation, the role of the in- 
tegral counter-term explicitly depends on the uv inte- 
gral cut-off of the integral (Q) in such a way that the 
three-bod y p arameter has a fixed value for all Q. In 
Refs. [6J, I63 a subtraction technique has been also in- 
troduced in order to regularize the atom-dimer scatter- 
ing problem. In these last references, the subtraction is 
made at zero momentum and therefore imposes the ex- 
act value of the atom-dimer scattering length. Thus, this 
regularization scheme does not correspond to the nodal 
condition of Eq. (|119l) . It is interesting to write down the 
equation for the atom-dimer scattering amplitude g{k) 
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obtained within this scheme: 

3g{k)a3 Soadas 



1 



^ ' 4 



(ka^y 
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1 + (kasY 



I^l/a3ik,u) 



2u^a 



2^2 



1 + u'^aj 



giu), (126) 



and to compare it with Eq. (J125I) . 



B. Heteronuclear trimers 

As a consequence of interesting predictions [1, ]1Q, IitI , 
l66l - [73J . heteronuclear systems play an important role in 
ultra-cold physics. In this section, the eigenequation for 
heteronuclear trimers is considered without external po- 
tential. However, it is worth noticing that an external 
potential can deeply change the physical properties of the 
system as it is the case for example in systems of mixed 
dimensions ^,[70]. Particles (labeled by i) can be ei- 
ther bosons or fermions and have a mass 771^ which can 
take two possible values: rrii = m or rrii = M . Bosons 
(respectively fermions) of mass M are denoted by B (re- 
spectively, by F) while bosons (respectively, fermions) of 
mass m are denoted by b (respectively, by /). Fermionic 
particles are supposed to have two possible internal states 
denoted f and \.. The configurations studied here are thus 
constructed from the set of particles {b, B, /^, /^, F^, F^}. 
The present study is also restricted to '2+1' few-body 
systems where the two-body interaction is non-vanishing 
for heterogeneous pairs of particles only. As an exam- 
ple in the case of (bF^) interacting pairs, other pairs like 
(bF^) or (6-B) are not interacting. For fermions neglect- 
ing (F-^F-^) or {f^fi) interactions is an exact assumption 
which follows from the Pauli principle. However, neglect- 
ing {BB) or [bb) interactions is not an exact hypothe- 
sis and only means that these interactions are negligible 
with respect to a heteronuclear interaction. Thus the 
true inter-atomic forces of the heteronuclear interacting 
pairs are in the vicinity of a s wave resonance while in- 
teraction for pairs of identical bosons is neglected. For 
convenience, the following notations are introduced for 
the different combinations of particles masses m.i : 



(127) 
(128) 

(129) 
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{ij) = 


rrii + rrij 


fJ-t 


(Jk) = 






j),m 


M(,j)M(^ki) 
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i\/r.. ., _L A/T,. .> 



'(u) 



'(kl) 



In this section, two identical atoms of mass M (particles 
1 and 2) interact with another one of mass m (particle 3). 
All the possible bound states can be thus deduced from 
the two configurations: (BBb) [or equivalently (BBf)] 
and (F^F^b) [or equivalently (F-^F-^/)]. An eigenequation 



M 




-f. 



FIG. 4: Schematic representation of the momentum coordi- 
nates used for the source amplitude associated with the con- 
tact condition between particles 1 and 3. 



for trimers is obtained for a negative energy E — E^, 



E. 



2fi 



<0, 



2,3 



(130) 



and if a dimer of binding energy £'2 exists then E^ < £'2 . 
In the center-of-mass frame, the source amplitude asso- 
ciated with the pair (23) is: 



231 |q2^3\ 



({rnis# 



(27r)^(5(ki+k2+k3)F(ki). (131) 



The other source amplitude with the momentum coor- 
dinates represented schematically in Fig. (|4|), is deduced 
from Eq. (|131|) by using the exchange symmetry between 
particle 1 and particle 2: 

({l"}|Sr') = si3(27r)^<5(ki + k2 + ka) F(k2), (132) 

where S13 = 1 when particles 1 and 2 are bosons and 
S13 = —1 if they are fermions. For simplifying the nota- 
tions, it is useful to introduce the mass ratio y defined 

by 



M 



y 



M2,3 



7713 AI + m 



(133) 



Assuming that the e ^- limit is well defined, the 
eigenequation is obtained from the contact condition for 
the pair (2 : M; 3 : m) in Eq. pUTj) where |^(o)) = 0: 






d^l 



(27r)'" u2 + fc2 + 2yk.u -|- q^ 



F(u) 



(134) 

In Eq. ()134p the identity d^k2d. — d^k2 has been used 
(the contact condition is performed at a fixed value of 
K23) and E^^ is the coUisional energy: 



El°' = £3 



2 1,2 



K'k 



2^1,(23) 



<0. 



(135) 



A detailed study of the three-body bound states in 2D 
has been performed in Ref. '67<], and in this subsection 
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FIG. 5: Critical mass ratio (■^), for the threshold of ap- 
pearance of an Efimov spectrum in each partial wave I of 

Eq. (Hazii. 



the discussion is centered on the 3D case. By using the 
rotational symmetry of the kernel in Eq. (|134p , one can 
fix an arbitrary direction: Cq and expand the source am- 
plitude F(k) in terms of partial waves as: 



F(k)=^Pz(ek-eq)Fz(fc). 



(136) 



1=0 



Each component -Fl(fc) verifies the integral equation 



~.«=^o, 



ynk 



2yku 



(137) 
where Qi is a Legendre function of the second kind. One 
can notice that the standard STM Eq. ()llip is obtained 
by setting formally S13 = 2 and y = 1/2 in Eq. (J137I) . 
For a sufficiently large value of the mass ratio, the Efi- 
mov effect appears in non zero partial waves [zl, UM ■ ^^ 
in the three-boson case, in this regime and in the limit 
of large momentum k Eq. (jl37p supports a pair of solu- 
tions Fi{k) - fc^2±<s (g g ^y The critical values of the 
mass ratio of the Efimov theshold in the I wave [denoted 
{^)i ] are plotted in Fig. ([5]). They have been deduced 
from Eq. (|137p in the Appendix [Al and the results found 
with this method coincide with the ones computed in 
Refs. i,[i3,[zl- 



Heteronuclear tetramers 



Recent theoretical [73|,|76|-|80| and experimental [8lll82| 
progress has been achieved in the four-body problem with 
ultra-cold atoms. In this section, a zero-range eigenequa- 
tion for four-particle bound states or "tetramers" is con- 
sidered. The binding energy is denoted E4 and 



2„2 



E4 = - 



2/7 



1,2 



(138) 



where E4 < E3 and/or E4 < 2E2 if a trimer and/or dimer 
exists. A particle (i) is characterized by a momentum k^ 
and a mass m^ equals to m or M. In the center-of-mass 
frame, the source amplitude associated with the pair (12) 
can be written as: 



({l''}|Sr') = (2vr)^<5(^ K)F{Ki2, ^34). (139) 

n=l 

The integral equation satisfied by the function F is ob- 
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FIG. 6: Schematic representation of the momentum coor- 
dinates used for the source amphtude associated with the 
contact condition between particles 1 and 2 for a system of 
fermions with two-mass components. 



tained from Eq. (|10ip (here (ij) — (1, 2)). A specific con- 
figuration is represented schematically in Fig. (j6]). The 
summation on the left hand side of the integral equa- 
tion is composed of source terms which are deduced from 
Eq. (|139l) by using the permutation symmetry and the 
statistics of the particles. In order to have a general 
equation for the different possible configurations, a sta- 
tistical factor denoted Snp is introduced for each pair of 
particles. Depending on the system, for two interacting 
particles n and p, the statistical factor is Snp = ±1 or 
Snp = for noninteracting particles. Using this notation, 
the different source amplitude are written in Tab. (JH. 
Different possible configurations of the two-mass compo- 
nent system and the corresponding statistical factors are 
listed in Table |lll Equation (jlOip can be simplified to 



d^ki 

(2^ 
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p>2 



F{Knp,kk 



'h'^F{u,v) 



n2 



2/ii,2rz3(£;£fV 

(140) 
In Eq. (dini) the 



where u = K12, v = k^^^. 

momentum k is related to the collisional energy by 

^cof^-^'«V(2Mi,2),thus: 
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^^^v'. (141) 
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Pair (ij) 


Permutations 


({r}isr> 


(12) 




4 

(27r)"5(^k„)F(A:i2,fc34) 


(13) 


2^3 


4 
Sl3(27r)^5(^k„)F(Ki3,fe24) 


(14) 


2^4 


4 
Si4(27r)"5(^k,0F(Kl4,fe32) 


(23) 


1 ^3 


4 

S23(27r)^5(^k„)F(A:23,fcl4) 


(24) 


1 o4 


4 
S24(27r)"5(^k,0F(A:24,fe3l) 


(34) 


lo3, 2o4 


4 
S34(27r)^5(^k,0F(A:34,fcl2) 



TABLE I: Source terms for each pair (ij) deduced by using 
the permutation symmetry of the four-body wave function 
from the source term of the interacting pair (12) in Eq. p39|) . 
The statistical factors Sij are given in Table HIl 



Configuration 


(mi, 7712,7713,7714) 


(Sl3, Sl4, S23, S24, S34) 


BBBB 


(m, m, 171,171) 


(1,1,1,1,1) 


if or b)BBB 


(m, M, M, M) 


(1,1,0,0,0) 


BbBb 


{M,m,M,m) 


(0,1,1,0,1) 


FthF^h 


{M,m,M,m) 


(0,-1,-1,0,1) 


BftBfr 


{M,m,M,m) 


(0,-1,1,0,-1) 


if or b)F^FrF^ 


{m,M,M,M) 


(-1,-1,0,0,0) 



TABLE IL Different possible 4-body configurations. The no- 
tation is as follows: B (respectively F) means that the atom 
is a boson (respectively fermion) of mass A/, 6 (respectively 
/) means that the atom is a boson (respectively a fermion) 
of mass 777. The fermions have two possible internal states f 
and J,. For each configuration, interaction is non vanishing 
only between one type of heterogeneous pair. The statistical 
factors Sij appears in the integral equation p42|) . 



Pair [ij) 


Vij 


(Tij 


(13) 


7774U 7774k 

T7 h V - — 

M3,4 A/2,4 


/'M1,2 

V '^1 


M3,4\ 
"73 / 


u V + k 


(14) 


7773U 7773k 

~T7 — 1""^+ Tr~ 

A/3,4 J1J2,3 


/'M1,2 


M3,4\ 
7774 / 


u + v + k 


(23) 


7774U 7774k 

h V -1 

Af3,4 Afi,4 


/M1,2 
\ 777i 


M3,4\ 
7774 / 


u 4- V — k 


(24) 


7773U 7773k 

T7 V-l- — 

Af3,4 Afl,3 


/'M1,2 
V rill 


M3,4\ 
1773 / 


u — V — k 



TABLE IIL Coordinates appearing in Eq. (|142 
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FIG. 7: Ground state branch for the three bosons of mass Af 
interacting with one impurity of mass tti in 2D as a function 
of the mass ratio. 



of 2D ground tetramers are computed numerically as a 
function of the mass ratio of the interacting particles in 
the BBBb configuration. Calculations are restricted to 
s wave tetramers by using the ansatz 



Using the notations defined in Table IIIII the integral 
equation which encapsulates the Bethe-Peierls contact 
condition can finally be written as: 



F(u,v) = F(u,i;,6') with 6* = Z(u, v). 



(143) 






iD 






F(k,vi3) 



2 






|-S24^7 



F(k,vi4) F(k,V23) 

K^ + 0-^4 K^ + 0-23 

F(k,V24) ^ ^ F(-u,k) 1 ^ -fi^F(u,v) 

(142) 

In the particular case where the system is composed of 
two heterogeneous dimers made of fermions {F-^ f^F-^ f^) , 
one recovers the same integral equation as in Eq. (16) 
of Ref. [iJI. In 3D depending on the statistics of the 
particle and of the mass ratio, Eq. (|142l) is not in gen- 
eral well defined as a consequence of the Efimov- Thomas 
effect. For systems composed of four particles in 2D, 
bound states of particles of same mass has been consid- 
ered in |77| . In the present work, the binding energies 



In order to check the numerical computation, the partic- 
ular case of four identical bosons which has been already 
obtained by several authors [7a, [Tg has been considered. 
In this configuration, two bound states have been found 
with binding energies which are close to already pub- 
lished results: E4/E2 — 197 and 24, to be compared with 
Ei/E2 = 197.3 and 25.5 in Ref. [t^ or with Ei/E2 = 194 
and 24 in Ref. [73| • Results for three identical bosons in- 
teracting resonantly with another particle is shown in 
Fig. ([7]). In 3D for the same bBBB configuration, simi- 
larly to the three-boson original STM equation, Eq. (|142p 
does not constitute a well-defined problem. The nature 
and properties of cigenstates of this system in the zero- 
range limit are yet unsolved. The bF^F^Ff or /jftfiJ^t 
configurations have been studied recently in Ref. |73| and 
a pure four-body Efimov effect (i.e., without a three-body 
Efimov state) involving a four-body parameter has been 
found for a mass ratio of 13.384 < M/m < 13.607. 
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VI. CONCLUSIONS 

In this paper, it has been shown how the zero-range 
approach can be formulated for different problems in a 
unified framework without technical intricacy for all di- 
mensions D < 3. In the few-body problem, general in- 
tegral eigenequations are obtained straightforwardly in 
the momentum representation. Efimov effect is deduced 
from the Skorniakov Ter-Martirosian equation by using 
a nodal condition and a subtracting technique which is 
close to regularizing schemes used in the Effective Field 
Theory. 

Various open issues remain to be solved for few-body 
systems in the limit of zero-range forces. For example, 
the question whether or not the 3-body parameter is 
enough for descri bing 4-bosons properties in 3D is still 
under debate [7^, ISOl • For heteronuclear few-body sys- 
tems, taking into account the external trapping (which 
is a natural concern in regards to experiments) makes 
the problem even richer. The external potential differs 
for different atomic species and this degree of freedom 
allows one to explore realistic intermediate situations be- 
tween the fully three-dimensional homogeneous configu- 
ration and the limit of mixed dimensions where inter- 
esting predictions have been already made in the three- 
body case [6^, [T^] ■ The present formalism can be also 
used to model systems where interaction occurs between 



three spin components like in ^Li experiments |83l - [86| . 
Another open issue concerns the possible occurrence of 
Efimov states for more than four particles in 3D. The 
four-body Efimov phenomenon has been found to occur 
for a specific interval of mass ratio in the 3+1 fermionic 
problem in Ref. |73| . Therefore, one can wonder whether 
or not a pure Efimov effect in a A'^ -|- 1 fermionic system 
can also occur and at which critical mass ratio for A^ > 3. 
Despite the difficulty of this problem, the integral equa- 
tion of the system in the limit of zero-range forces can be 
easily deduced from the present formalism and is given 
below as a concluding remark. Each of the A^ polarized 
fermions of mass M labeled by i (1 < « < A^) interact 
with only one impurity (particle i ^ N + 1) of mass m. 
In the center-of-mass frame, the source amplitude for the 
pair A^, (A^ -I- 1) can be written as: 
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N.,{N+i)\\(^N^N+l 



}isr 



xF(ki,k2...kw-i), (144) 



and other source amplitudes [S^^'^^) {i = 1- ■ ■ N —1) 
are deduced from this ansatz by using the fermionic 
statistics. The general integral equation for this problem 
which encapsulates the Bethe-Peierls asymptotic condi- 
tion is then obtained from Eq. (I101[) : 
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i^(ki,k2,- ••kAr_2,k 
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2m \l^i=l "■' 
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tAT-l) 



TD{Ecoi + iO+) 
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d^k 
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(27r 



i-D 



*("'(ki,---kAr), (145) 



2M 



E£'<^? 



(£ 



Af-1, 



where ^coi - ^ ^ 2M ^i=i "-* - 2(A/+m) y^i=i 
^(") is a complementary solution associated with the free 
Hamiltonian at energy E in the center-of-mass frame (it 
is equal to zero for E < 0) and F is antisymmetric under 
the exchange of two coordinates (k^ , kj ) . 
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Appendix A: THRESHOLDS FOR THE 
'l+2'-BODY EFIMOV EFFECT 



In this appendix, the values of the mass ratio for which 
an Efimov effect occurs are deduced from Eq. (|137p . To 
this end, zero energy solutions at unitarity jaaj =00 are 
considered (solutions for finite energy have the same high 
momentum behavior). In this regime, Eq. (I137P is scale 
invariant and this allows one to search for power law solu- 
tions: Fi{k) — fc'^+'~^, where t/ is a function of the mass 
ratio. For convenience, the following dimensionless pa- 
rameter t is introduced: 



t = arcsin(j/) = arcsin 



M 



M- 



(Al) 



Thus, the scattering amplitude is /3(Kk) = l/(A:cost) 
and Eq. (|137p gives for each partial wave I an eigenvalue 
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equation: X{v,l,t) = where, 



\[v,l,t) ~ cost — 



TTsint 



dzz 



v+l-l 



Qi 



l + z^ 



2zsmt ^ 
(A2) 
The integrand in Eq. (IA2[) is positive, hence for a Bose 
(respectively Fermi) statistics only even (respectively 
odd) values of I can support a solution. For x > 1, 
Qi{x) = iP,(x)ln(f±i) + Wi{x) where Wi{x) is a poly- 
nomial of order I. The polynomial Wi does not contribute 
in the integral of Eq. (jA2|) and \{v,l,t) can be thus ex- 
pressed in terms of a sum of functions of the form: 



lil.t) 



dz z^ In 



z^ -I- 2z sin t + 1 
z^ — 2zsint + 1 



(A3) 



with < t < f and |3fi(7)| < 1. In Eq. (jMl), -^(7,*) can 
be calculated explicitly as: 



I{l,t) 



27rsin(7i) 
7C0s(^)- 



(A4) 



In the / wave, Efimov states are characterized by a power 
law such that v = —I -\- is, where s is real. For example, 
for two identical bosons interacting with another particle 
(si3 = 1), in the s wave sector. 



X{i',l = 0,t) = cost 



sm{i't) 



'sin(t)cos(^) 



(A5) 



Equation (|A5|) admits solution of the form v = is for all 
values of the mass ratio and s -> for i — >■ 0. For higher 
I waves, Efimov states appear above a critical value of 
the mass ratio. The threshold is obtained by searching 
the value of the parameter t = i?"^'* such that: 



and the threshold is obtained from the equation 
i™* =tan(if'*)- f sin2(i5;"*), which gives the critical 

mass ratfo {^Yi% - 13.6 found in Refs. d [zl. For 
higher partial waves, the critical values of the mass ra- 
tio ( — ), found by this method coincide with the ones 
computed in Ref. [17[ where the hyperspherical method 
was used. Results are reported in Fig. ^. 

In the regime where an Efimov effect occurs, the regu- 
larizing technique can be achieved by fixing the asymp- 
totic phase shift which is a function of the three-body 
parameter (denoted kI) in the partial wave /: 



Fi{k) ex 7:7 sin 

fc-i-oo K 




(A8) 



A set of solutions satisfying Eq. (jASp can be filtered from 
Eq. p37p by imposing the nodal condition: 



mi 



'D 



where k^^f 






,pen. (A9) 



Exact zero-range theory is obtained in the limit where 
p -> cx). Universal results are recovered if the node k^J'f is 
chosen at a large value as compared to all the low energy 
scales {kl'f » l/|a|,(7). Using the subtracting scheme of 
the three-boson problem, Eq. (IA9|) can be also incorpo- 
rated into the integral equation Eq. (|137p : 



lim X (1^,1,^1"^) =0. 



(A6) 



In the p-wave sector which concerns the case of two 
identical fermions interacting with another particle 
(si3 = — 1) one gets 



\{i',l ^ l,t) = cost - 



(y + 2) s\n{vt) - v s\n[{v + 2)t] 
2v{v + 2) sin^ t cos tcosi^) 



iATL 



I /-oo 




(AlO) 
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